Topological Order, Dimerization, and Spinon Deconfinement 
in Frustrated Spin Ladders 
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We consider topological order and dimer order in several frustrated spin ladder models, which are 
related to higher dimensional models of current interest; we also address the occurrence of fraction- 
alized phases with deconfined spinon excitations in these models. Combining results obtained with 
both analytic and numerical methods, we discuss how the occurrence of dimerized or fractionalized 
phases are dictated by the system's geometry. 
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I. INTRODUCTION 

Frustrated spin systems give rise to a wealth of inter- 
esting behavior; hence, they have attracted considerable 
attention. 1 Typically, the low- lying excitations in spin 
systems are magnons. However, there has been particu- 
lar interest in identifying fractionalized phases, where the 
magnon "breaks apart" into more fundamental objects. 
When doped, such phases would give rise to metallic 
states which do not fall within the Fermi liquid paradigm 
— the elementary excitations in these phases would not 
have the quantum numbers of an electron. However, 
finding systems which exhibit fractionalized phases has 
proven to be extremely challenging. 

Motivated largely by Anderson's original resonat- 
ing valence bond (RVB) suggestions for the high-T c 
cuprate superconductors *r substantial effort has focussed 
on searching for fractionalized phases in models having 
short-range RVB ground states 4 A related approach has 
been to consider dimer models, where only short-ranged 
valence bonds with specified dynamics are considered.— 
While lacking long-range order in the conventional sense, 
i.e., lacking a local order parameter, these short-range 
valence-bond ground states have a subtle form of order, 
namely topological order. 5 It is now appreciated that 
a precise characterization of fractionalized phases is via 
its topological order.— However, whether a fractionalized 
phase actually occurs depends strongly on the system's 
geometry. More specifically, the system's geometry must 
allow for "liquidity" in the spectrum of states^ Many sys- 
tems with short-range RVB ground states do not exhibit 
such liquidity — they would prefer to dimerize, rather 
than exhibit a fractionalized phase. 8 

In this work, we address the occurrence of fraction- 
alized and dimerized phases in several frustrated spin 
ladder models. Ladder models provide a unique test- 
ing ground, as powerful analytic and numerical tech- 
niques from one-dimensional physics can be utilized. In- 
deed, ladders models have allowed controlled calcula- 
tions to investigate topological order ; 9 ' 10 : 11 the occur- 
rence of dimer order ^ and the occurrence of fractional- 
ized excitations^ in spin models. Furthermore, ladders 



models have allowed for controlled calculations demon- 
strating that pairing and, in particular, d x 2_ y 2 pairs 
could arise when these spin models are dopedJ^ The 
models we consider are related to higher-dimensional 
systems of current interest; investigations of the one- 
dimensional analogs are particularly relevant, in light 
of recent work which showed that the mechanism giving 
rise to fractionalized excitations is the same in both one 
and two dimensions J£ However, besides being a testing 
ground and illustrating the types of possible behaviors, 
these ladder models are interesting in their own right, as 
there are a number of materials that are well described 
by ladder models J£ 

In large regions of parameter space, the models we con- 
sider have short-range RVB ground states. More specif- 
ically, in large regions of parameter space these models 
have ground states which are continuously related to the 
ground states of the so-called rung-singlet phase or the 
Haldane phase. [Typical configurations in these ground 
states are shown schematically in Figs. [2] and Q] in Sec. II.] 
While these phases have nearly identical properties, their 
ground states differ in a subtle way, namely in their topo- 
logical order. With frustrating interactions, these states 
could be tuned to become degenerate, and then one has 
the necessary liquidity for the deconfinement of spinon 
excitations. However, as will be seen below, depending on 
the model's spatial symmetries, this liquidity and spinon 
deconfinement could be preempted by dimerization. 

The rest of the paper is organized as follows. In Sec. II, 
we describe the ladder models considered in this work — 
the cross-coupled, zigzag, and diagonal ladders — and 
we recall some of their known properties. In Sec. Ill, the 
exact ground states of the models are discussed (along 
certain lines in parameter space). Sec. IV contains a dis- 
cussion of topological order in these models and some 
comments on the low-energy excitations. In Sec. V, 
we analyze the models in the limit of weak interchain 
couplings, using bosonization and renormalization group 
techniques. In Sec. VI, we present numerical results 
obtained via the density-matrix renormalization-group 
(DMRG) algorithm; we discuss these numerical results 
in light of the results from the previous sections. Finally, 
in Sec. VII we summarize and present some concluding 
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remarks. 



II. THE MODELS 

We begin with two antiferromagnetic spin-1/2 Heisen- 
berg chains, each described by the Hamiltonian 



H 
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where Sj^ is the spin operator at site I on chain i. We will 
be interested in coupling the chains together in various 
ways, such that the resulting models have different spa- 
tial symmetries. In this work, we will consider only anti- 
ferromagnetic couplings. We start by coupling the chains 
together so that the resulting models are not frustrated. 
We then include frustrating interactions and investigate 
their influence on the properties of the models. 

The simplest way to couple the chains together is by 



so that the full Hamiltonian is 



(2) 



= E4 
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which is depicted in Fig. [TJ We will refer to this as an 
ordinary ladder. This model has received considerable 
attention. Indeed, it can be thought of as a strip of a 
two-dimensional square lattice; hence, its properties have 
been investigated to give clues as to the physics occurring 
in the high-T c cuprate superconductors^ 
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FIG. 1: The ordinary ladder. 

When J± is large (Jj_ ^> Jy), the ground state is essen- 
tially a product of rung singlets with a gap to the excited 
states — the gap is due to the energy necessary to break 
a singlet bond. When J± ~ Ju, it has been shown that 
the energy gap persists; the ground state is well described 
by a short-range valence-bond state, a typical configura- 
tion of which is shown in Fig. It has been established 
that the entire region < J± < oo is, in fact, contin- 
uously relatedi 18 i 19 As this entire region of parameter 
space is related to the regime J± Ju where the ground 
state is a product of rung singlets, it is often referred to 
as the rung-singlet phase. Incidentally, the dominance of 
rung-singlet bonds can be measured by determining their 
weight in the ground state p s = (I/AT) J2h=i(^\ s i s l\^) 



where |\&) is the ground-state wave function and S l — 

(1/V2) [| T) ; (1) 1 1) ; (2) - 1 | T)| 2) ] • A state belongs to the 
rung-singlet phase if p s > 1/4. 20 
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FIG. 2: A typical valence-bond configuration in the rung- 
singlet ground state of the two-leg ordinary Heisenberg ladder. 

We also consider coupling the spins to their next- 
nearest neighbors on the opposite leg of the ladder. This 
coupling is described by 



F X = £J X (S«<> + SP>.S« 



i 

so that the full Hamiltonian is 

2 



H C = Y,W + H 



(4) 



(5) 



i=l 



the resulting model is shown in Fig. [3] 
Ju 




FIG. 3: Ladder with diagonal couplings. 

This system is know n 21 1 22 to have a gapped spectrum 
of spin-1 magnons. Furthermore, the entire regime < 
Jx < oo is continuously related to the Haldane phase 2 ^ of 
the spin-1 chain. This can be understood by considering 
the point Jx = J\\ — here, the low-energy spectrum of 
(© is equivalent to that of the 5 = 1 Heisenberg spin 
chain] 24 ' 25 To see this, we start with the Hamiltonian of 
a spin-1 chain H = Ju Si ■ Sj+i, which is known to 
be in the Haldane phase. When the spin-1 operator on 
site I is represented as a sum of two spin-1/2 operators, 



Si = Sj +Sj , one obtains |(SJ) with Jx = Ju. Since the 
total spin of each rung commutes with the Hamiltonian, 
the eigenstates can be classified by the total spins on the 
rungs. It has been shown that in the low-energy part 
of the spectrum all rungs are in their triplet (5 = 1) 
state, and hence the same Haldane gap appears in the 
ladder model as well. In this representation the Haldane 
state of the spin-1 chain can be described rather well by 
short-ranged valence bonds between neighboring rungs, 
a typical configuration of which is shown in Fig. 2] 

Finally, we consider the ladder model shown in 
Fig. HJa) , where the interchain coupling is given by 



(2) 
21 + 1 



(6) 
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FIG. 4: A typical valence-bond configuration in the ground 
state of the ladder model shown in Fig. [3] 

Previous work has established 25 that this model is in the 
same universality class as the model in ([3]). Hence, re- 
moving half of the cross couplings in an appropriate way 
from Eq. |5} does not change the universality class. Part 
of our motivation for considering this model is because 
it is a representation of the ll N p = 2 diagonal ladder" 
shown in Fig. [5f b) ; it is a minimal model to study di- 
agonal stripes, which have been observed in the high-T c 
cuprate superconducting material I^-^Sr^CuC^^ and 
the nickel oxides La2Ni04.i29^ and Lai-ajSr^NiO/t^ 
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FIG. 5: (a) Ladder with a diagonal coupling between every 
second spin, (b) N p = 2 diagonal ladder. 

As mentioned above, we are interested in the effect of 
frustration on these models. In what follows, we analyze 
three models obtained by introducing frustrating inter- 
actions to the models described above. (Again, we will 
consider only antiferromagnetic interchain couplings in 
this work.) Our motivation for doing so is because the 
resulting models have different spatial symmetries. As 
discussed below, these symmetries play a crucial role in 
determining the properties of the models. 

1. Cross-Coupled Ladder 

The first model we consider has both H± and Hx 
present simultaneously. The resulting model with Hamil- 
tonian 

2 
i=l 

is shown in Fig. [SJ For the rest of this work, we will 
refer to this model as a cross-coupled ladder. Notice 
that this spin model is invariant under translation by 
a single site; it is also invariant if the legs of the lad- 
der are interchanged: leg 1 <-> leg 2. Part of our mo- 
tivation for the study of this model is its relationship 



to the two-dimensional model with nearest-neighbor and 
next-nearest-neighbor exchange couplings, often referred 
to as the J\ — J2 model — (J7J) is a one-dimensional strip 
of this model. Furthermore, {7J can be thought of as a 
chain of edge-sharing tetrahedra. 30 As there are a number 
of materials described by spin models on corner-sharing 
tctrahedra, i.e., pyrochlores, 31 it is not unreasonable that 
materials with edge-sharing tetrahedra can be realized. 




FIG. 6: The cross-coupled ladder. 



The model in ([7]) has been investigated in a num- 
ber of works, but its phase diagram and properties are 
still under debate. When J± is the dominant interchain 
coupling, the system is continuously related to the or- 
dinary ladder (with Jx = 0); hence, the system is in 
the rung-singlet phase. On the other hand, when Jx is 
the dominant interchain coupling, the system is continu- 
ously related to Eq. §5§ , and the system is in the Haldane 
phase. For weak interchain couplings, previous analytic 
treatments 1 ^ suggested a first-order transition between 
the rung-singlet and Haldane phases when J± ~ 2Jx- 
Numerical results on the mode l 11 ! 32 were consistent with 
a first-order transition for both weak and strong inter- 
chain couplings. Recent numerical work, however, has 
suggested the transition is actually continuous when the 
interchain coupling is weak, becoming first-order only 
for stronger interchain couplingi 20 ' 33 Furthermore, re- 
cent analytic worki£ has argued there is a spontaneously 
dimerized phase in between the rung-singlet and Haldane 
phases, rather than a direct transition between the two 
phases as observed in previous works. 

2. Zigzag Ladder 

The second ladder model we consider contains again 
both rung and diagonal couplings, but only half the di- 
agonal couplings of the cross-coupled ladder are present, 
as shown in Fig. [7J This model is often referred to as a 
zigzag ladder] its Hamiltonian is 

2 

H z = Y / H i o ) +H ± +H Z , (8) 

»=i 

where 

ff z = ]Tj 2 S^-S^. ( 9 ) 

I 

This spin model is invariant under translation by a single 
site. However, unlike the cross-coupled ladder which is 



4 



invariant under the interchange of the legs, this model 
lacks that symmetry. Part of the motivation for con- 
sidering this model comes from the two-dimensional tri- 
angular lattice — © is a one-dimensional strip of the 
triangular lattice. The spin- 1/2 Heisenberg model on a 
triangular lattice has been of considerable recent inter- 
est, motivated largely in part by the discovery of the 
triangular lattice material CS2CUCI4 and, in particular, 
to evidence that this material exhibits a two-dimensional 
fractionalized phase with dcconfined spinons. 34 However, 
besides being a toy model for understanding spin sys- 
tems on a triangular lattice, this model is relevant to the 
quasi-one-dimensional material SrCu0 2 . 35 




FIG. 7: The zigzag ladder model. 



The properties of the zigzag ladder are well known 
when J±_ or J 2 vanishes, and also along the line J±_ = J 2 . 
When J± or J 2 vanishes, the model reduces to the ordi- 
nary ladder. When J± = J 2 , the zigzag ladder is equiv- 
alent to the spin-1/2 frustrated Heisenberg chain with 
nearest-neighbor coupling J± and next-nearest neighbor 
coupling J11 . This chain model is known to be critical for 
Ji > J Xc (J±c = ^11/0.241)^^2. being in the same uni- 
versality class as the antiferromagnetic spin-1/2 Heisen- 
berg chain. The spinons of the spin-1/2 chain acquire a 
gap (but they remain deconfined), and the ground state 
becomes doubly degenerate for J± — J 2 < J± c ~ At the 
Majumdar-Ghosh (MG) point 39 J± = J 2 = 2J\\, the two 
degenerate ground states have a simple form in the ther- 
modynamic limit, each consisting of decoupled singlets, 
as shown in Fig. [8j It has been shown that these two 
ground states can be continuously related to the rung- 
singlet and Haldane phases4£ For even weaker interchain 
couplings J± — J 2 < 2Jy, incommensurate oscillations 
appear in the short-range correlationsi 41 i 42 i 43 i 44 i 45 



3. Diagonal Ladder 

Finally, we consider the model obtained when the cou- 
plings described both by ([2]) and ([6]) are present simulta- 
neously. The model is shown in Fig. Ufa) and is described 
by the Hamiltonian 

2 

H D = J2H + H ± + H d . (10) 

i=l 

For the rest of this work, we will refer to this ladder 
model as a diagonal ladder. Unlike the previous two 
ladder models which were invariant under translation by 
a single site, (|10j) is invariant under translation by two 
sites. Furthermore, similar to the zigzag ladder, (|10|) is 
not invariant under the interchange of the legs. However, 
it is invariant under the combined operation of leg inter- 
change followed by a translation by a single site. Interest- 
ingly, this ladder model can, in fact, be transformed into 
one which is invariant under translation by a single site, 
but lacking inversion symmetry. This is accomplished by 
interchanging the sites on every second rung of the lad- 
der, S 2/ <-> S 2/ ; the resulting model is shown in Fig. [9] 
We are unaware of work which investigated the role of 
frustrating interactions in this model. 




FIG. 9: (a) The diagonal ladder model, (b) Alternative rep- 
resentation of the diagonal ladder, obtained by interchanging 
sites on every second rung. 



III. EXACT GROUND STATES 




FIG. 8: The two degenerate ground states of the zigzag ladder 
at the Majumdar-Ghosh point. 



An interesting feature of the models considered in this 
work is that their exact ground states have a simple 
form for a certain regime of parameters. The key to es- 
tablishing this is decomposing the system into a set of 
"triangles" , 4 - Here, we review the argument and discuss 
the ground states that arise. 

Consider three s = 1/2 spins arranged on a triangle. 
Suppose the spins on sites 1 and 2 are coupled by an ex- 
change coupling Ji; these two spins are coupled to the 
spin on site 3 by an exchange coupling J. The Hamilto- 
nian for the triangle is 



-^triangle = Jl&l ' S 2 + J (Si + S 2 ) • S 3 . (11) 
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We are interested in the case where both J and Ji are 
antifcrromagnetic; hence, the ground state is a doublet. 
There are two ways of achieving this: (a) a singlet across 
sites 1 and 2 with a free spin-1/2 on site 3, or (b) a 
triplet across sites 1 and 2 added to the spin-1/2 on site 
3 to form a spin-1/2. Let denote a singlet between 
sites i and j. Then, using that 



(S, + S J )[i,j]=0, 



(12) 



the energy of the state in (a) is found to be E a — — 3 Ji/4; 
the energy of state (b) is = — J + Ji/4. We see that 
the state in (a) — the state with a singlet across sites 1 
and 2 — is the lowest-energy state if J < J\ . 

Using information from the previous paragraph, we 
now address the ground states of the ladder models con- 
sidered in this work. We begin by considering the cross- 
coupled ladder model. In this model, four triangles can 
be assigned to each rung; the Hamiltonian can be written 
as 



H, 



cc 



i 



H, 



(i) 
i+ 



H, 



(i) 



H, 



(2) 
1 + 



H 



(2h 



(13) 



with 



We now consider the zigzag ladder. It can be decom- 
posed into two triangles per rung: 



I 



(4 2) 



H 



«1 



where 



H, 



(2) 
1 + 



(2) 



H, 



(i) _ i 



i- 



- 2 J ^ S l 



(!) . S (2) 



/„s (2) -s (2) 



(18) 



1 7 C(l) q(2) 

2 J 2°l ' Oj+i ) 

(19) 

1 7 c(2) c(l) 
2 J 2»; ' S i-1 • 



When Jii = J2/2, the terms in (|19p can be written in the 
form of l|ll[): the energy of the state in (jX5j) saturates 
the lower bound in (fT7|) for Jj_/2 > J||. Hence, (fTS]) is 
the exact ground state for J|| = J2/2 with J±/2 > J\\^ 
Note that another exact ground state of the zigzag ladder 
follows by symmetry. When J|| = Jj_/2 and J2/2 > 
J|| the ground state is again a product of singlets, but 
this time they are formed diagonally between neighboring 
rungs 



40 



(20) 



H, 



(i) 
1+ 



\J±S l • Sj 



(2) 



+ 2 J II S ; 



(i) 0(1) 



j (2) 



rr(l) _ 1 7 c(l) c(2) , 1 7 c(l) q(l) , 1 T C 



1+1 1- 2 ,J ^l 
(!) , 1 7„a(2) 



(2) 



1 7 q(1) 



(1) «(2) 
/ 

s(2) 



\ J \\ S l 



(2) e (2) 



+ 1 J II S Z 



(2) 



j(2) 



(1) 



(1) 
l+l ' 

(1) 
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(2) 



J (2) 
'(-1 • 

(14) 



We now go on and consider the diagonal ladder. The 
diagonal ladder can be decomposed into two triangles 
per rung, but the triangles have different orientations for 
even and odd rungs: 



H n = 



H 



(2) 
(20- 



-H, 



(i) 

(2/+1) 



-H, 



(i) \ 

(21+1)-) 



(21) 



When Jii = Jx, the terms in 
form of IJTTJI- Then, using (TT 



14)) can be written in the 
, it follows that the state 



w=n[iM2] 



(15) 



where the spins on each rung form a singlet, is an exact 
eigenstate of the Hamiltonian in ([13]) with energy 



E 



:J±N. 



(16) 



We would like to determine if and when (|15|) is the ground 
state. To do so, consider each of the triangles in (fT4|) 
individually. If J±/2 > Jn, the lowest-energy state of 
each triangle (when considered individually) has a singlet 
between the spins on the same rung. It is known that if 
a Hamiltonian can be written as the sum of terms, the 
ground-state energy cannot be smaller than the sum of 
the lowest energies of its constituents. For the ground 
state of (Tl^|) . this gives the inequality 



E >--J ± N, 



(17) 



where now 
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H, 
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(21+1)- 
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+" 2 J 2» 2 ; 
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21 
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j(2) 



J 2I 



JnS 



(2) e (2) 



\°2l 



J 2l-1 



(1) «(2) 



2( + l 
, 1 j c(l) 

+ 2 J 2» 2 ;+i 

- I T ^ 

— gJJ_o 2 , +1 • 

, 1 j o(2) 

+ 2 J 2» 2 ;+i 
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21+1 
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When Jii = J2/2, the terms in (|22|) can be written in the 
form of HJ); and thus (fTS"]) is the exact ground state for 
J|| = J 2 /2 with J ± /2 > J||. 



IV. TOPOLOGICAL ORDER AND SPINONS IN 
LADDER MODELS 



where N is the number of rungs. Hence, we see that (TTJ 
saturates the bound in (fTTj) when Ji/2 > Jn; (fTS]) is the 
ground state for Jn = Jx with Ji/2 > Ji ] 30 i 47 



^From the discussion in Sees. II and III, the ladder 
models considered in this work have similar ground-state 
properties in extended regions of parameter space - 
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their ground states are described by a collection of short- 
ranged valence bonds, separated by a gap to the excited 
states. However, the ground states of these ladder mod- 
els, in fact, differ in a subtle way, namely in their topo- 
logical order. More specifically, the number of valence 
bonds crossing an arbitrary vertical line is always even 
in the rung-singlet phase, while the number is always 
odd in the Haldane phase. This can be seen explicitly 
in the configurations shown in Figs. [2l [4] and [8] Hence 
a topological number Q can be defined by the parity of 
the number of short-range valence bonds crossing an ar- 
bitrary vertical line^iS This Q, which is either even or 
odd, is a good quantum number for short-range valence- 
bond states since the Hamiltonian has finite matrix ele- 
ments only between configurations with the same Q. For 
long-range valence-bond states, however, the Hamilto- 



nian mixes the Q = even and Q = odd configurations; 
hence, no such topological distinction is possible. 

It is worth noting that for open boundary conditions 
Q = odd ground states have spin-l/2's localized at the 
ends of the ladder, while Q = even states do not. As 
can be seen from Figs. 0]and[Hb), these end spins oc- 
cur for topological reasons. They are analogous to the 
edge states in the quantum Hall effect; in general, the 
presence of such edge excitations is a signal of nontrivial 
topological order.— 

It has also been pointed out^S that the topological or- 
der of the valence bonds is related to the "hidden order" 
present in two-leg Heisenberg spin ladders, namely string 
order, analogous to the string order in antiferromagnetic 
spin-1 chains.— This string order is detected by the two 
string order parameter o 9 ' 10 i 50 



Oeve„ = - ( (S&u + S&) exp (for E (Sf +1 , 1 + S£ 2 )J (S? +M +S£ 3 )). (23) 



l=i+l 



It was showni* that a slightly modified version of these 
string order parameters can be written in terms of the 
order and disorder fields of two Ising models — Oodd 
can be written in terms of the order fields <7\ and 172; 
Ocvon can be written in terms of the disorder fields \i\ 
and fj,2- Hence, O dd and Oovcn cannot be nonzero si- 
multaneously, as they are dual to each other. It was 
observed that models with Q = odd ground states have 
Oodd 7^ 0, while models with Q = even ground states 
were observed to have O e ven 7^ 0i lo i n Hence, the string 
order parameters detect the topological order of the va- 
lence bonds. 

This topological order has important consequences for 
the excitation spectrum. In these models with short- 
range valence-bond ground states, the simplest excitation 
is generated by breaking one of the valence bonds, pro- 
moting it to a triplet. An interesting and important ques- 
tion is how this excited state propagates. More specif- 
ically, does the triplet propagate coherently, or does it 
"break apart" so that the individual spins forming the 
triplet — referred to as spinons — propagate indepen- 
dently? In both the rung-singlet and Haldane phases, 
it is known that the triplet propagates coherently. The 
reason for this can be understood looking at Fig. [Poland 
counting the number of valence bonds crossing a verti- 
cal line. When a valence bond is broken in a state with 
Q = even, if this triplet breaks apart it leaves a string of 
valence bonds with the "wrong" topology in the interme- 
diate region. This gives rise to an increase in the local 



energy, which is proportional to the distance between the 
two spinons; as a result, the spinons are confined into a 
(gapped) spin-1 magnon. (A similar situation occurs in 
the Q = odd Haldane phase.) 




FIG. 10: Spinons separating topologically distinct regions. 

/From Fig. [lOl one also sees that if the two topologi- 
cally distinct ground states can be made degenerate, one 
can expect spinons to be deconfined. This is because 
the string of "wrong" valence bonds between the two 
spinons would not give an increase in energy. Thus, the 
degeneracy provides the "liquidity" necessary for spinon 
deconfinementii Indeed, as discussed in Sec. II, this is 
known to happen in the zigzag ladder due to the frus- 
trating interaction. It is reasonable to expect similar 
phenomena could occur generically — frustrating inter- 
actions could tune the topologically distinct rung-singlet 
and Haldane ground states to be degenerate, so that 
spinons are deconfined and propagate as elementary ex- 
citations. While this expectation is reasonable, as will be 
discussed in detail below, there are other possible phases 
- namely dimerized phases — which could intervene. 
Moreover, we will see that the system's geometry plays 
a crucial role in determining whether the system may 
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dimerize or not. 



V. WEAK-COUPLING ANALYSIS 

It is useful to consider the physics in the limit where 
the interchain coupling is weak, as controlled analytic 
calculations are possible. More specifically, we start with 
two decoupled spin- 1/2 chains and consider the inter- 
chain coupling as a perturbation. Provided that there 
is no phase transition, the low-energy Hamiltonian de- 
duced in the weak-coupling limit is valid even at strong 
coupling, albeit with renormalized parameters. This ap- 
proach has been utilized in various other works; besides 
providing a detailed understanding of the properties of 
the two-leg ladder and its relation to the spin-1 chain^ 
it has been effective at uncovering and illucidating the 
phenomena that can arise in these systems i 12 i 13 ' 51 i 52 ' 53 i 54 

We begin this section by briefly describing the formal- 
ism, mainly to establish our conventions. (More detailed 
accounts can be found, e.g., in Ref. 1 551.1 We go on to de- 
duce the effective low-energy Hamiltonians for the mod- 
els considered in this work, and then discuss the physics 
contained in these effective Hamiltonians. In particular, 
we discuss the role of the various irrelevant operators 
that arise. These operators were ignored in most previ- 
ous works, but were recently argued to affect the physics 
qualitatively in some situations^ As will be elaborated 
on below, quantum fluctuations can give rise to several 
different behaviors. To understand the physics of these 
quantum fluctuations, we use the renormalization group 
(RG) and examine the behavior of the system under a 
change of scale. 



A. Formalism 

The low-energy properties of the spin-1/2 chain are 
described by an SU(2)i WZW model with Hamiltonian 



H = — J dx(3 R -3 R + 3 L -3 L ), (24) 

where the velocity v is related to Jm, and 3 R and 3l 
are currents satisfying the SU(2)i Kac-Moody operator 
product expansion (OPE) 



: J%{z) :: jg(«, 



: J£(z) :: jf (w) 



1C 



-Jl{w). 



2(z~w) 2 ' 2{z-w) RK J ' 

(25) 



a/37 



2{z-w) 2 2{z-w) 



JIM: 



with z = vt + ix and z = vt - ix. The SU(2)i WZW 
model has a single primary (matrix) field g(z,z) of di- 
mension (1/4,1/4). (g is the field appearing in the a- 
model representation.) Physical operators of the spin- 
1/2 chain are given by combinations of the components 



g — the staggered magnetization n a (z, z) and the dimer- 
ization e(z, z): 

n a (z,z) = Ti[a a g(z,z)} , e(z,z) = Tr[ 5 (z,z)] , (26) 

where the {<r Q } are Pauli matrices. The n a and e fields 
have leading short-distance behavior 



n a {z) :: n^(w) : = 



: e(z) :: e(w) 



\z — w\ ' 
1 



\z — w\ ' 



(27) 



: n a (z) :: e{w) : = 0. 



Furthermore, being linear combinations of the compo- 
nents of g, their OPE's with the currents are 



: J&z) :: n^w) : = 

: J£(z) :: n^w) : = 

: J%{z) :: e(w) : = 

: J£(z) :: e(w) : = 



2(z 


— w) 




i 


2(z 


— Tu) 




i 


2(z 


— w) 




-i 


2(z 


— To) 



n a (w). 



(28) 



[In Eqs. d2Tj) and (J28|), n a (z) = n a (z,z) and e(z) = 
e(z,z).] 

At low energies, the spin operator on chain i, Sj f \ can 
be written as 



j(') 



1 



M 



a 2tt 



(3 iR (x) + 3 iL (x)) + (-iy——n i (x) , (29) 

27ra 1 ' z 



where a is an ultraviolet regulator (a ~ the lattice spac- 
ing), and M is a nonuniversal 0(1) constant. In what 
follows, we will need to know the behavior of the fields 
upon translation by a single site. ^From the invariance 
of the spin operator S, under translation we deduce 
how 3i R (x), 3iL{x), and rii{x) transform. Then, as all of 
the terms on the right-hand-side of (|28p must transform 
in the same way upon translation, we deduce how €i(x) 
transforms. Hence, we arrive at 

3iR/L(n) -> 3m/l(x), (30) 
n-i(x) — > -nj(x) , 6i(x) — > ~£i(x) . 



B. Operators and Phases: Topological Order, 
Dimer Order, and Deconfined Spinons 

As described above, starting from two spin-1/2 chains, 
we are interested in the fate of the system upon turning 
on an interchain coupling. In this subsection, we discuss 
the various operators that can arise in the low-energy 
Hamiltonian(s); we discuss the physics that these opera- 
tors give rise to. 
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The low-energy Hamiltonians are dictated by symme- 
try. The models we are considering are all SU(2) sym- 
metric; hence, all operators appearing must transform as 
SU(2) scalars. Therefore, if only relevant and marginal 
operators are considered, the operators 

ei and e 2 , (.\t 2 , ni • n 2 , nid x n 2 , eid x e 2 , 
(Jin + Ji L ) ■ n 2 , (J2H + Jax) • n x , (31) 

(JlR ' J2L + J2R ' Jli) i (Jlfl ' JlL + J 2 L ' Jail) 

have the potential of appearing in the low-energy Hamil- 
tonians. Furthermore, we will not consider the terms 
Jir • J2A and In, ■ J 2 l- Even though these terms have 
dimension two (and, hence, are marginal), they couple 
excitations moving in the same direction on both chains 
and give rise only to small quantitative corrections. As 
discussed below, other symmetries of the various models 
will further restrict the operators allowed in their low- 
energy Hamiltonians. 

The physics is determined by the operator which flows 
to strong coupling first under the RG. When the ni • n 2 
term determines the physics, the resulting phases have 
gapped magnons with spinons being confined^ [The 
ni • ri2 term, in fact, gives rise to a strong confining po- 
tential which binds together spinons from the two chains 
to form a gapped magnon.] Furthermore, these phases 
have topological order; information about this topolog- 
ical order is contained in the sign of the ni • n 2 term's 
coefficient. Another possibility is if the physics is deter- 
mined by the (Jir ■ 3 2 l + J2R ■ Jil) term. When this 
occurs, it has been shown that the resulting phase is a 
fractionalized phase with deconfined spinonsJ^ Finally, if 
the physics is determined by the e\ and e 2 or €162 term(s), 
the result is a spontaneously dimerized phase. 

There are other operators in Eq. (|3"Tj) . However, these 
operators do not determine the phases that arise, though 
they modify the properties within a particular phase. In 
particular, (Jir ■ Jil + J2L ■ ^2r) is the marginally ir- 
relevant operator present in the two spin- 1/2 chains.— 
Furthermore, ni • d x n 2 and e\d x e 2 are believed to give 
rise to incommensurate correlations ! 13 i 54 Finally, as will 
be seen below, the (Ji^ + Jil) • n 2 and (3 2 r + 3 2 l) ■ ni 
are always subleading to ni • n 2 and 6162 and, hence, do 
not determine the phases that arise. 



C. Effective Low-Energy Hamiltonians 

Here we deduce the effective low-energy Hamiltonians 
for the frustrated ladder models considered in this work. 
We then derive renormalization group (RG) equations 
describing how the parameters in the low-energy Hamil- 
tonian evolve under a change of scale. Finally, we derive 
the initial values of the parameters from the microscopic 
models in the limit of weak interchain coupling. 



1. Cross-Coupled Ladder- 
Resides being SU(2)-symmetric, the cross-coupled lad- 
der is invariant under translation by one site and also 
inversion about the center of the ladder: leg 1 <-> leg 2. 
Translation by one site forbids ei and e 2 from appear- 
ing alone; it also forbids the operators (Ji_r + Jil) • n 2 
and (3 2 r + 3 2 l) ■ ni. Inversion about the center forces 
operators from leg 1 and leg 2 to appear symmetrically; 
in particular, it forbids the operators nid x n 2 and e\d x e 2 
from appearing. These symmetries constrain the effective 
low-energy Hamiltonian to have the form 

H = J cfe [7 (3 1R ■ 3 1L + 3 2L ■ 3 2R ) (32) 

9 P 
+A (3 1R ■ 3 2L + Jil • 3 2R ) + - ni • n 2 + - e\t 2 

a a 

To understand the physics contained in (|32[) , we use the 
RG and investigate the behavior under a transformation 
of scale. Using the OPE's in ||57J, and (EH}, we 

deduce the RG equations for the parameters to be^ 



n 

dX 

H 

dg 

dl 

d/3 



9- ]pg + ^g- 7; A/3, 
P + \iP-\\g. 



(33) 



^From the structure of the RG equations, we see that if 
either g or (3 is nonzero, both ni • r\ 2 and £162 will be 
generated upon renormalization. However, if both g = 
and (3 = 0, these operators will not appear. Hence, one 
would expect the relevant ni • xi 2 or e\t 2 terms to deter- 
mine the physics under most situations. If both of these 
terms are suppressed, the physics would be determined 
by the marginally relevant (Jlr • J2L + Jil • J2_r)- 

To deduce the values of the parameters in (|32| . we 
insert {29]) into {2J and (gj); we obtain 



f dx r. T „ T N M 2 
J j2^i {J± ~ 2Jx) — ni ' n2 



„ T M- v-^ a""" 
- 2J X Yl 7^T n i^ "2 

m=l v ; 

+ (J x + 2J X )(3 1R ■ 3 2L + 3 1L ■ 3 2R ) 



(34) 



A lattice Hamiltonian gives rise to a continuum field the- 
ory with an infinite number of operators.— Here, besides 
ni-ri2 and (Jij? ;• J2L + Jil • 3 2 r), an infinite number of ir- 
relevant operators of the form ni • d x m n 2 appear. ^From 
Eq. (|34|) . the relevant operator ni • n 2 determines the 
physics for generic values of J\ 1 and J± . When the ■ n 2 
term is suppressed, it appears the physics is determined 
by the (Jir ■ 3 2 l + Jil ■ J2_r) term. However, things are 
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more subtle — the irrelevant ni • (9 2m n 2 can generate the 
relevant ri! • n 2 and £162 when integrated outJ^ Indeed, 
focussing on the leading irrelevant operator ni • d 2 n 2 , the 
term 



2Jx) 



d 2 x\ 



(J 



li? ' •>2L 



Jli 



>2R) 



x J x M 2 a I — X2 ni • <9 2 n 2 



(2tt) 2 



(35) 



appears in the partition function at second order. Using 
the OPE's in fjgl), 



Jx(Jl + 2 J x )Af 2 
(27r) 2 a 



1 3 
-ni • n 2 - -e x e 2 



(36) 



is generated in the low-energy Hamiltonian. Hence, to 
leading order we deduce 



A = 



J± + 2 Jx 
(2tt) 2 



p = 



3 ( J± + 2J X ) J X M 2 
8tt (2tt) 2 



.9 = 



(Jl - 2 J X )M 2 J X M 2 Ji. + 2 J x 



(27T) 2 



4tt 



(37) 



deduce the RG equations for the parameters to be 



~dl 
~dl 

dg 

dl 
dp 

~dl 
dgi 
dl 

<m_ 
di 



7 "2<? +431-4/31, 
A 2 + .g 2 - .g/3 + ig 2 - -gift. , 
.g- ^.g + Ag- -A/3, 
/3+^7/3-^Ag, 
-^73i + A51 - -A/3i , 

2 7/?1 "2 A31 - 



(39) 



As with the cross-coupled ladder, ni • n 2 and eie 2 will 
be generated upon renormalization unless both g = 
and /3 = 0. Similarly, ni ■ dn 2 and e\de 2 will appear 
unless both g\ = and /3i = 0. Hence, as with the 
cross-coupled ladder, one would expect the ni -n 2 or eie 2 
terms to determine the physics under most situations; if 
they are suppressed, the physics would be determined by 
(Jii?. • J2L + Jil • J2fl)- 

To deduce the values of the parameters in Eq. (f3"5|) , we 
insert Eq. (J25J) into Eqs. © and ©; we obtain 



As mentioned above, the relevant operator ni-n 2 deter- 
mines the physics for generic values of J|| and Jj_. Now 
when g — > in (|37[) . the relevant eie 2 is present (from 
integrating out the irrelevant operators); one may ex- 
pect it to determine the physics. However, its coefficient 
is significantly smaller than the (Jir ■ J 2 l + Jil ■ J2_r) 
term. Therefore, a subtle competition between the two 
interactions should be expected. This competition will 
be investigated in detail in our numerical calculations. 



2. Zigzag Ladder 

Like the cross-coupled ladder, the zigzag ladder is in- 
variant under translation by a single site. However, it 
lacks the symmetry of inversion about the center of the 
ladder. Hence, the operators from leg 1 and leg 2 do 
not have to appear symmetrically — now the operators 
nid x r\2 and ei<9 x e 2 are allowed. The effective low-energy 
Hamiltonian has the form 

H = J dx 7 (3 ir ■ Jil + J2L ■ J2L>) 

+ A (J lfl • J 2L + 3il ■ 3 2R ) (38) 

+- ni • n 2 + - eie 2 +gi n r d x n 2 + Pi eid x t 2 ■ 
a a J 

As with the cross-coupled ladder, we investigate the 
physics of Eq. (gHJ) using the RG. Using the OPE's, we 



H ± +H z = 
M 2 — 



+ 



dx 
(2^ 

-(- 

m! V2 



{{Ji 



M 2 

J 2 ) ri! • n 2 

a 



(Ji + Ji) (Jli? ■ JlL + J 2 L ' J 2 i?)} 



(40) 



As we saw with the cross-coupled ladder, the lattice 
Hamiltonian gives rise to an infinite number of operators 
in the field theory. Similar to the cross-coupled ladder, 
the relevant ni • n 2 determines the physics for generic 
values of J|| and J±. Furthermore, when the ni ■ n 2 
term is suppressed, the irrelevant operators may gener- 
ate terms which determine the physics. Interestingly, for 
this model all the terms of the form ni • <9 2m n 2 have their 
coefficient proportional to (J± — J 2 ). Hence, all the op- 
erators of the form ni • (9 2 ™ l n 2 are fine-tuned away along 
the line J± = J 2 ; only terms of the form ni • 9 2m+1 n 2 are 
present. As there are only derivatives of odd power, they 
are unable to generate terms like ni • n 2 or eie 2 . There- 
fore, for this model the irrelevant operators do not change 
the physics (though they give rise to small quantitative 
changes, e.g., in the size of energy gaps); to leading order 
we deduce the parameters in Eq. (f3"B")) to be 



(Jj_ - J 2 )M 



J± + J 2 
(2tt) 2 ' a ~ (2tt) 2 
(J x + J 2 )M 2 



9i = 



/3 = 0, 
01=0. (41) 
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3. Diagonal Ladder 



insert (l29l) into flOl); we obtain 



Like the previous two models, the diagonal ladder has 
translation invariance. However, unlike the cross-coupled 
and zigzag ladders which are invariant under translation 
by a single site, the diagonal ladder is invariant under 
translation by two sites. Furthermore, while this model 
is not invariant under inversion about the center of the 
ladder (leg 1 <-> leg 2), it is invariant under inversion com- 
pounded by translation by one site. Due to the invari- 
ance under translation by two sites, the operators t\, £2, 
(Jifl+Jli) - n.2 and (J2ii+J2L)'Hi are now allowed. How- 
ever, the symmetry of inversion compounded by transla- 
tion by one site constrains them to appear in the combi- 
nation (ei-e 2 ) and [{3\ R + J 1L ) ■ n 2 - (J 2 fl + 3 iL ) ■ n x ]. 
Hence, the effective low-energy Hamiltonian has the form 



H 



dx 



a 



(ei - e 2 ) + 7 (Jir • Jil + J 2 l ■ J 



2R) 



g (3 

+A (Jifl • J 2 l + Jil • J2i?) + - ni • n 2 + - eie 2 

a a 

H — 777 [(Jlr + Jil) • n 2 - (J 2 i? + J2L) ■ ni] 



(42) 



As before, we use the OPE's to deduce the RG equa- 
tions for the parameters; we find 



d- t 

~dl 
(IX 

~dl 

dg 

dl 

dp 

~dl 

da 

~dl 

dn 

~dl 



2 1 2 1 o2 1 2 

7 -2 9 + 2 + 2 a ' 
\ 2 +g 2 -gf3 1 

1 1 „ 

9- ^19 + *9- 2 X P~ K ' 
3 3 

-a+ -pa + pa, 

11 1 1 

-ac — — 7« + — A/t — <?k + -/3k; . 



(43) 



(2^{ (JW2) 



M 



(2 

M 2 00 _2m 

a ^ (2m)l 

(J± + J2) (Jlfl ' JlL + J 2 L ' J 2 fl) 



-ni • n 2 



(44) 



+ ^fy^ [(Jib. + Jil) • n 2 - (J 2 r + ^2l) • n x ] 



00 „2m 



m— 1 



(J 2 B+J2L)-a 2m ni ] } 



Interestingly, even though the term (£1 — 62) is allowed, its 
coefficient is zero. However, the relevant [(Ji_r + Jil) • 
ri2 — (J2_r + J2L) • ni] is always present; as mentioned 
above, it generates both the nin 2 and £i£2 terms. Hence, 
for this model the irrelevant operators give rise only to 
small quantitative corrections; the physics that can oc- 
cur is already contained in the relevant operators already 
present. Therefore, to leading order we deduce the values 
of the parameters in (|4"2"|) to be 



A = 



J± + J 2 



9 = 



(Jl - J 2 )M 2 



{2nf 



{3 = 0, 



J 2 M 



a = 0. 



(45) 



As mentioned above, the zigzag ladder is expected to 
have a fractionalizcd phase when g — > ; one can expect a 
subtle competition between £162 and (Jij?-J2L+JiL , J2i?) 
in the cross-coupled ladder. For this ladder model, as Ei£2 
is generated by the relevant [(Jir + Jil) ■ 112 — [3 2 r + 
^2l) • ni] term, it is reasonable to expect a regime of 
parameter space where £162 flows to strong coupling first 
and the system dimerizes. 



VI. NUMERICAL RESULTS AND DISCUSSION 

In this section, we present our numerical results for the 
phase diagrams of the models introduced in Sec. II; we 
discuss them in light of the analytical results presented 
in the previous sections. 



For this ladder model, the operator (£1 — £2) is allowed; 
if its initial coefficient is zero, this term will never be 
generated. Furthermore, we must have a = 0, g = 0, 
[3 = 0, and k = in order for all the relevant operators 
to be banished. In particular, if any of these are nonzero, 
the ni • n 2 and £i£2 terms will be generated. 

To deduce the values of the parameters in Eq. (|4"2"|) , we 



A. Numerical Method 

The numerical calculations have been performed on fi- 
nite ladders with open boundary condition (OBC) using 
the DMRG algorithm 59 with the dynamic block-state se- 
lection (DBSS) approach i 60 ' 61 We have set the threshold 
value of the quantum information loss x to 10 -8 and the 
minimum number of block states M m { n to 64. All relevant 
eigenstates have been targeted independently using four 
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to six DMRG sweeps until the entropy sum rule has been 
satisfied. The accuracy of the Davidson diagonalization 
routine has been set to 10~ 7 . 

Recently, it has been shown that quantum phase tran- 
sitions (QPTs) can be conveniently studied by calculating 
some measure of entanglement.2S.S2iSlSSiSS.SLSm2a.ZL2a 
In particular, the von Neumann entropy of a block con- 
taining a finite number of neighboring sites often gives a 
clear indication of a QPT, as anomalies appear in these 
quantities at the transition: the entropy exhibits a jump 
at a first-order transition or develops a cusp (with in- 
creasing N) at a continuous transition. It should be 
noted, however, that depending on how a first-order tran- 
sition is realized, one might have difficulty distinguishing 
it from a continuous transition. More specifically, if the 
two levels corresponding to the different ground states 
are already orthogonal in a finite-sized system, the en- 
tropy of a block will exhibit a jump when the two levels 
cross (in a finite-sized system). However, if the two lev- 
els are not orthogonal in a finite-sized system, the wave 
function, energy, and consequently the entropy of a block 
will vary continuously in any finite-sized system. In this 
case, the level crossing develops only asymptotically, and 
the jump in the entropy appears only in the N — > oo 
limit. 

In this work, we consider (i) s/, the entropy of the I th 
rung, (ii) sy+i, the two-rung entropy of the neighboring 

I th and (I + l) st rungs, (iii) the two-site entropy of 

the spins Sj and Sj^Vj on chain i, and also (iv) sat(Z), 
the entropy of a block formed by the left / rungs of a lad- 
der with N rungs. To avoid end effects, we compute si, 
sij+i, and s\ fl \ 1 in the middle of the ladder, for / = N/2 
or I = N/2 +1. As discussed above, one of our pri- 
mary interests is to identify dimerized phases which may 
arise and the concomitant breaking of translational sym- 
metry. The appearance of a columnar dimerized phase 
can be detected by considering the difference of two-rung 
entropies 

D a = 81+1,1+2-81,1+1 l = N/2, (46) 

Alternatively, taking the block entropy sn(1) of the left 
I rungs, one can consider 



D s = s(l) - s(l + 1) / = N/2 . 



(47) 



which tells how the block entropy of the left half changes 
when an extra rung is added. The appearance of stag- 
gered dimerization can be detected by considering the 
difference of two-site entropies on the two chains 



P - (s w -s (L > )-(s {2 > -s {2 > ) I -N/2 

(48) 

Further information about the phases that arise can 
be deduced by studying the length dependence of 
sn (l) i 58 ' 73 i 74 ' 75 i 76 For noncritical, gapped models, this 
quantity saturates to a finite value when I is far from 



the boundaries, while for critical systems 



sn(1) = g In 



~2N 




sin 




7T 





(49) 



where c is the central charge . 58 ' 77 Moreover, if the sys- 
tem's ground state is spatially inhomogeneous, oscilla- 
tions appear in sat(/)4 5 - Hence, the Fourier spectrum 



~s(k) 



1 N 

-Y 

N f-^ 



e~ lkl s N {l) 



(50) 



carries information 45 about the spatial inhomogeneity: if 
the amplitude of a peak at a nonzero wave number k* 
remains finite in the thermodynamic limit, this indicates 
a periodic spatial modulation of the ground state with 
wavelength A = 2tt /k* ', 

In what follows, we will often need to know the large- N 
behavior of various quantities. For any quantity A, the 
finite-size-scaling Ansatz 



A(N) = Ao + a/N? 



(51) 



has been used, where Aq, a, are free parameters deter- 
mined by a least-squares fitting procedure. 



B. Results and Discussion 

1. Zigzag ladder 

We first present results for the zigzag ladder. Since 
the phase diagram and properties of this model are well 
known, it serves as a test case for the other models. The 
Jx i J2 parameter space was explored by calculating var- 
ious entropy functions for J± and J2 satisfying 



J 1 



Ji = CJ n 



(52) 



for several values of C . In the numerical calculations the 
energy scale was set by taking Jm to be unity. 

We first show results obtained for C = 2. As seen in 
Fig. [Ill s i m the middle of the ladder is discontinuous at 
Jx = J2 = J\\ j indicating a first-order transition. Besides 
the discontinuity, the inset of Fig. [TTI shows a minimum in 
the entropy at J2 — 1.36 indicating a possible second 
transition. Moreover, there is also a minimum in s(k) at 
an incommensurate k* in the region 0.5 < J2/ J\\ < 1-36. 
However, on physical grounds, there is no reason to ex- 
pect a second transition in the zigzag ladder since any 
deviation from the J± = J2 line drives the system into 
the rung-singlet or Haldane phase. Indeed, this mini- 
mum is a finite-size effect — it originates from the two 
end spins of the Haldane phase. To substantiate this, 
notice that there is no sign of this minimum in the rung- 
singlet phase, while we know that the phase boundaries 
must be symmetric under the interchange of J± and J2. 
Furthermore, by attaching spin-l/2's to the ends of the 
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ladder with a strong anitferromagnetic coupling, the end 
spins can be eliminated^ When the calculations are re- 
peated with these extra spin-l/2's attached to the ends, 
the minimum in s(k) at an incommensurate k* disap- 
pears. 



2x32 
2x64 
1.3 -0- 2x128 




0.94 0.96 0.98 1 1.02 1.04 1.06 1.08 1.1 
J 2 



FIG. 11: (Color online) Single-rung entropy as a function of 
J-2 for Jj_ + J2 = 2J|| with J|| = 1 for three different lengths 
of the ladder. 




FIG. 12: (Color online) Single-rung and block entropies for 
J± + J2 — 10 J|| , as a function of J2 for ladders with 32, 64 
and 128 rungs. The insets show the behavior near J2 = 2 Jy , 
where the entropies vanish. 

Performing the calculations for couplings satisfying 
(152")) with other values of C, a behavior similar to Fig.fTTI 
was observed whenever the line J± = J2 was crossed at 
J± < J||/0.241. That indicates that the system under- 
goes a first-order transition (along the line J± = J2) 
for J± < J||/0.241. Different behaviors in the en- 
tropies were observed, however, for stronger interchain 
couplings. Fig. [T2l shows si (top panel) and s;,z+i (bot- 
tom panel) in the middle of the ladder for J±+ J2 = 10 J11 . 
In the rung-singlet phase si is small, as the spins on a 



rung are predominantly in a singlet state; in the Haldane- 
like phase si is close to In 4 ~ 1.386, since the valence 
bonds are formed predominantly between neighboring 
rungs. Notice, however, that si no longer has a jump 
at J± = J2 = CJ\\/2 (as it did in Fig. [TTjk rather, si and 
s/j + i have a discontinuity in their slope. When sn(1) 
was computed at the point J± = J2, it was found that it 
could be fit well with the form given in (|4T)|) with c = 1. 
These results are consistent with a continuous transition 
between the rung-singlet and Haldanc-like phases in this 
regime. They are consistent with the fact that, if the 
zigzag ladder is written as a single chain with nearest- 
and next-nearest-neighbor couplings, the chain is criti- 
cal in this regime with the low-energy physics being de- 
scribed by (J2U). 

The insets of Fig. [12] show si and s^+i in a region 
about J2 = 2 Jy ; we see that the entropy functions vanish 
at J2 = 2 Jy for any length of the ladder. This is because 
the exact ground state is a product of rung singlets [see 
Eq. (fl~5|) ] for J2 = 2J||. In fact, the entropies vanishes 
along the entire line J2 = 2Jy when J± > 2Jii, as (fl5)) 
is the exact ground state. Furthermore, computing the 
entropy of the two spins coupled along the diagonal by 



J2, S ; and Sj^-p this entropy was found to vanish for 
J± = 2Jy when J 2 > 2J||. This is because (|2T))) is the 
exact ground state along this line. More generally, due 
to the symmetry of the model, one obtains the same re- 
sults presented in Figs . [TT1 and [T2l if J± and J2 are inter- 
changed; instead of considering rung entropies, one must 
consider the entropies of diagonally coupled spins. 
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FIG. 13: (Color online) Phase diagram of the zigzag lad- 
der. The first- and second-order transition between the rung- 
singlet and Haldane phases are denoted by solid and dashed 
lines, respectively. The ground state along the dashed-dotted 
fine is an exact rung-singlet or diagonal-singlet valence-bond 
state. The dotted lines indicate the parameter values wher 
the DMRG calculation were done. 

Putting together these results, we arrive at the phase 
diagram in Fig. 1131 The transition line at J± = J2 is 



of first order in the weak-coupling limit, which becomes 
of second order at strong couplings. The entire region 
above (below) the line J± — J2 is continuously related to 
the exactly solvable line J2 = 2Jn (J± — 2Jy) and, hence, 
is in the rung-singlet (Haldane) phase. In terms of the 
effective low-energy Hamiltonian discussed in the previ- 
ous section, the physics in these phases is determined by 
the ni • 112 term; in the rung-singlet (Haldane) phase, 
its coefficient is positive (negative). As discussed in 
Sec. IV, the rung-singlet (Haldane) phase has Q — even 
(Q = odd) topological order, and the elementary excita- 
tions are gapped magnons with spinons being confined. 
Along the J± = J 2 line, the two topologically distinct 
ground states become degenerate, spinons are deconfincd, 
being domain walls between the two topologically dis- 
tinct (and energetically degenerate) ground states. As 
mentioned in the previous section, at weak coupling the 
physics is determined by the (Ji_r • J2L + ^2R ■ Jil) term 
(with the ni -d x n2 and t\d x ti terms giving rise to incom- 
mensuration). Looking at Eq. (TiD|) . this fractionalized 
phase occurs because the geometry of the zigzag ladder 
fine tunes away an infinite number of operators which 
could cause dimerization. 



2. Cross-Coupled Ladder 

As with the zigzag ladder, we calculated s; and sy+i 
for J± and Jx satisfying (f52|) with various values of C. 
Fig. [14] shows results obtained for intermediate values of 
the interchain couplings, C = 2 and C = 3. Here, we see 
the entropy functions display a finite jump, indicating 
a first-order transition. The inset shows that s; vanishes 
(independent of the ladder's length) at Jx = J\\ ■ This oc- 
curs because (fT5j) is the exact ground state for Jx = J\\ 
with J± > 2J||. It was shown in Ref. |32| that (fl"5|) is, in 
fact, the ground state for Jx = J\\ and J± > 1.401 J11. 
Our numerical results are in agreement with this predic- 
tion. 

Besides the finite jump at Jx = 0.8, Fig. [TJ] also shows 
that the two-site entropy possesses a minimum around 
Jx ~ 1; furthermore, weak incommensurate oscillations 
were found to appear in sn(1)- Similar to the zigzag lad- 
der, these are finite-size effects due to the end spins in the 
Haldane phase. Indeed, when calculations are repeated 
for a system in which spin-l/2's are attached to the ends 
of the ladder (to freeze the end spins), the block entropy 
Sn(1) saturates for shorter chains and the minimum in 
s(k) at an incommensurate k* disappears. The ampli- 
tude of the remaining negative peak in s(k) at k* = n 
was found to vanish in the thermodynamic limit. Hence, 
no spatial inhomogeneity develops along the transition 
line. 

The numerical results show different behavior when 
C < 1.3. si as a function of Jx is shown in Fig. [Pol for 
C = 1. The ground-state wave function is continuous 
in the weak-coupling regime. As longer and longer lad- 
ders are considered, si exhibits a sharper and sharper 
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FIG. 14: (Color online) Two-rung entropy of the cross- 
coupled ladder as a function of Jx for Jj_ = 2J« — Jx for 
various system sizes. The dashed line is a guide to the eyes. 
The inset shows the same but for the J± = 3J\\ — Jx- 




0.9' 

0.2 0.4 0.6 0.8 1 

J x 

FIG. 15: (Color online) One-rung entropy of the cross-coupled 
ladder model for Jx + Jx = J11 , for several system sizes. The 
solid line is a polynomial fit. The inset shows the finite-size 
scaling of the position of the maximum of the entropy. 

maximum (bounded from above by In 4) , developing into 
a cusp at Jxc/J\\ = 0.355(3). Such behavior is sugges- 
tive of a continuous transition. Recent numerical works 
reported a continuous transition in the weak-coupling 
regime ) 20 i 33 in disagreement with analytic results^ as 
well as previous numerical calculations iii However, the 
analytic results are expected to be reliable in the weak- 
coupling regime. As was checked by the DBSS procedure, 
a reliable extrapolation of the gap requires calculations 
on longer ladders and keeping a significantly larger num- 
ber of block states than was available in Ref. |20. The 
same holds for the entropy. Moreover, as was shown 
above, spurious effects can arise due to end spins. When 
the calculations were repeated by attaching spin-l/2's to 
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the ends of the ladder; si behaves somewhat differently. 
As shown in Fig. 1161 instead of an abrupt change in slope, 
a jump now seems to develop; its position scales to the 
same value of Jx c obtained in Fig. 1151 This behavior sug- 
gests that the transition is of first-order. Based on these 
considerations, we believe the transition is, in fact, first- 
order. Although the singlet ground-state wave function 
is continuous in the weak-coupling regime for finite-sized 
systems, a crossing with the next singlet level (which we 
found to be at relatively high energy for weak couplings, 
far from the exactly solvable line) may develop in the 
N—>oo limit, and the asymmetric cusp in Fig. [15] may 
develop into a jump. 




0.33 0.34 0.35 0.36 0.37 0.38 0.39 0.4 0.41 

J.. 



FIG. 16: (Color online) Same as Fig. [15] but with spin-1/2 
particles attached to the two ends of the ladder with strong 
antiferromagnetic couplings. The inset shows the finite-size 
scaling of the position of the jump in the entropy. 

A further interesting feature of the model for weak in- 
terchain coupling is shown in Fig. |T7j More specifically, 
the two-rung entropy s^+x measured for I — N/2 ex- 
hibits two well separated peaks at Jx Cl (N) and Jx C2 {N) , 
but it exhibits a single peak for I = N/2 + 1. The differ- 
ence D S (N) is also finite in the region between Jx Cl and 
Jxc 2 - This could suggests the existence of a columnar 
dimer pha se in a narrow range of couplings, as predicted 
in Ref. 12. However, the two peaks merge in the N — > oo 
limit and the width of the putative dimer phase shrinks 
to zero. As shown in the inset of Fig.[T71 a finite-size scal- 
ing analysis gives the same critical value as the one-rung 
entropy. The same behavior was observed in the calcula- 
tion along the line Jx = 0.2Jh confirming the findings of 
Ref. HI 

We have also computed sjy(0 and its Fourier transform 
s(fc). \s(k)\ was found to have an extra peak at k* = 7r, 
besides the one at k = 0. However, \s(k* = n)\ was 
found to vanish in the large- N limit, indicating that the 
ground state is always spatially homogeneous. Similar 
behavior was found for other values of J± and Jx in the 
weak-coupling regime. Furthermore, we computed the 
staggered dimerization [see Eq. (|48p] and we found that 




FIG. 17: (Color online) Two-rung entropy of the cross- 
coupled ladder model for J± + Jx = Jn for various system 
sizes. The solid line is a polynomial fit. The inset shows the 
finite-size scaling of the positions of the peaks of the entropy. 



it also vanishes in the N — » limit. This provides strong 
evidence that, at least for weak interchain coupling, there 
is no intermediate columnar or staggered dimer phase 
between the rung-singlet and Haldane phases. 



O 1 st order 
- Exactly solvable 




FIG. 18: (Color online) Phase diagram of the cross-coupled 
ladder. The symbol x denotes the transition point calculated 
in Ref. |33|. A denotes the point where the ladder model is 
equivalent to a spin-1 chain. The ground state is an exact 
rung-singlet state along the dashed-dotted line. 

Putting together these results, we obtain the phase 
diagram for the cross-coupled ladder, shown in Fig. 1181 
The entire region above the transition line is continu- 
ously related to the exactly solvable line Jx/J\\ = 1 and, 
hence, is in the rung-singlet phase. The entire region be- 
low the transition line is continuously related the point 
Jx — 1, Jj_ = and, hence, is in the Haldane phase. 
As with the zigzag ladder, the physics in these phases is 
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determined by the ni • ri2 term in the low-energy Hamil- 
tonian — these phases have topological order, and con- 
fined spinons. As discussed in the previous section, at 
the transition there is a subtle competition between the 
(Jir ■ J2L + Jil • ^2R) and £162 terms, which give rise to 
a fractionalized and dimerized phase, respectively. How- 
ever, the entire transition line appears to be of first-order, 
with no evidence for an intermediate dimerized phase be- 
ing found. Hence, our results suggest that a dimerized 
phase does not appear in this model. In constructing the 
phase diagram, we made use of a duality relationship of 
the model. Similar to what was described for the diag- 
onal ladder in Sec. II, one can interchange the spins on 
every second rung: S^/- When this is done for this 

ladder model, another cross-coupled ladder is obtained 
but with Jy and Jx interchanged. This implies that en- 
ergies and, in particular, energy gaps satisfy 

E(J h J x , Jx) = E(J x ,J x , J,,). (53) 

Scaling by J|| , one obtains 

E(J x /J\\,Jx/J\\) = (Jx/J\\) E(J x /Jx,J\\/Jx). (54) 

3. Diagonal Ladder 

As with the other models, we computed si and s;,;+i for 
J_i_ and J2 satisfying (|5"2")1 with various values of C. Here, 
in contrast to the other two models, we find a dimerized 
phase intervening between the rung-singlet and Haldane 
phases. 




J 2 

FIG. 19: (Color online) Two-rung entropy of the diagonal 
ladder model for J± + J2 = J\\ for various system sizes. The 
solid line is a polynomial fit. The inset shows the finite-size 
scaling of the positions of the peaks of the entropy. 

Fig. [19] shows sy + i for Jj_ + J 2 = Ju. While the 
behavior appears similar to Fig. I17( the N — ► 00 limit 
is drastically different. Indeed, while a single transi- 
tion was found for the cross-coupled ladder, the inset of 



Fig.[T9lshows the two peaks do not collapse in this model. 
Hence, Fig. Q1J] suggests the system undergoes two dis- 
tinct transitions at J2 Cl /J\\ — 0.459 and Jic 2 l J\\ = 0.563, 
with a columnar dimerized phase between the two transi- 
tion lines. To substantiate this, we have computed sn(1) 
and subsequently \s(k* = 7r)|, the entropy difference be- 
tween neighboring plaquettes D s (shown in Fig. I20|) . and 
also the energy difference between neighboring plaque- 
ttes. All of these quantities were found to scale to a 
finite value in the region between J2 Cl and J2 C2 . Similar 
behavior was found for other values of the parameters in 
the weak-coupling limit. Hence, we conclude a columnar 
dimerized phase does, in fact, exist between the rung- 
singlet and Haldane phases in this model. 




FIG. 20: (Color online) The dimerization entropy of the diag- 
onal ladder model for J± + J2 ~ J\\ for various system sizes. 
The solid line is a polynomial fit. The left panel shows the 
finite-size scaling of the D s for various J2 values. 

When the calculations were repeated for Jx+ J% > 2J|| , 
a single first-order transition was obtained. It is worth 
noting that due to the rather small values of the gap, the 
asymptotic behavior can be seen only for long ladders. 
As before, convergence can be accelerated in the Haldane 
phase by attaching spin-l/2's to the ends of the ladder 
(to pin the end spins). 

Fig. 1211 summarizes our finding for the phase diagram 
of the diagonal ladder. The entire region above the tran- 
sition line is continuously related to the exactly solv- 
able line J2I — 2 and, hence, is in the rung-singlet 
phase. The entire region below the transition line is in 
the Haldane phase. As with the previous two ladder 
models, these phases have topological order with con- 
fined spinons; the physics in these phases is determined 
by the ni • 112 term in the low-energy Hamiltonian. How- 
ever, contrary to what was found for the cross-coupled 
and zigzag ladders, from the discussion in the previous 
section, it is reasonable to expect a regime in between 
the rung-singlet and Haldane phases where the diagonal 
ladder dimerizes (i.e., where the physics is determined by 
the ei£2 term). This is, indeed, found to be the case — a 
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FIG. 21: (Color online) Phase diagram of the diagonal ladder. 
The first and second-order phase transition points are indi- 
cated by the circle and square symbols, respectively. The ex- 
act rung-singlet ground state is indicated by a dashed-dotted 
line. 



narrow, but extended region is found where the ground 
state is dimerized. The first order transition is replaced 
in the weak coupling regime by two second-order transi- 
tions. 



VII. CONCLUDING REMARKS 

In this work we considered several frustrated spin lad- 
der models, which are related to higher-dimensional mod- 
els of current interest. In large regions of parameter 
space, these models have short-range RVB ground states 
with topological order; they are adiabatically related to 
the ground states of the rung-singlet or Haldane phase. 
We investigated the role of frustrating interactions on the 
models, addressing in particular how the transition be- 
tween phases with different topological order occurs. In 
the simplest direct transition takes place along the 

line where the even- and odd-topology phases become de- 
generate. While the elementary excitations of the topo- 
logically ordered phases are gapped magnons, spinons be- 
come deconfined along the transition line and a fraction- 
alized phase is obtained. Alternatively, the transition 



may occur in two steps, with an intermediate phase hav- 
ing broken translational symmetry. In this case, spinons 
always remain confined. 

An important observation from our analysis is the 
strong "desire" for broken-symmetry phases to arise at 
the transition between the rung-singlet and Haldane 
phases. Indeed, we saw that spin models typically give 
rise to an infinite number of operators that could cause 
dimerizationJ^ These operators arose in the diagonal lad- 
der, and a dimerized phase was seen to appear in the 
phase diagram. These operators also arose in the cross- 
coupled ladder; nevertheless, in agreement with Ref. l33l . 
we found no evidence for dimerized phases in our nu- 
merics. The RG equations suggest that this occurs due 
to a subtle interplay/competition of quantum fluctua- 
tions. Although no dimerized phase appears in the cross- 
coupled ladder studied in this paper, in this delicate sit- 
uation even small perturbations are likely to drive the 
cross-coupled model into dimerized phases as shown in 
Ref. [79|. The zigzag ladder was an exception — the 
model's geometry fine-tunes away the infinite number of 
operators which could cause dimerization. 

We believe these results give an outlook into the 
physics of higher-dimensional systems. In particular, our 
results show the importance of a system's geometry in 
achieving the necessary liquidity for fractionalized excita- 
tions to occur. Hence, our results illustrate why fraction- 
alized phases are hard to come by — fractionalized phases 
are delicate objects, requiring some level of fine-tuning. 
Even on lattices where the necessary fine-tuning occurs 
(such as the triangular lattice), small perturbations due 
to, e.g., spin-phonon coupling^ or ring exchanges are 
likely to drive the system into a dimerized phase. 
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